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Abstract 

We develop the technique of the variational approximation for solitons in two directions. First, 
one may have a physical model which does not admit the usual Lagrangian representation, as 
some terms can be discarded for various reasons. For instance, the second-harmonic-generation 
(SHG) model considered here, which includes the Kerr nonlinearity, lacks the usual Lagrangian 
representation if one ignores the Kerr nonlinearity of the second harmonic, as compared to that 
of the fundamental. However, we show that, with a natural modification, one may still apply the 
variational approximation (VA) to those seemingly flawed systems as efficiently as it applies to 
their fully Lagrangian counterparts. We call such models, that do not admit the usual Lagrangian 
representation, semi-Lagrangian systems. Second, we show that, upon adding an infinitesimal tail 
that does not vanish at infinity, to a usual soliton ansatz, one can obtain an analytical criterion 
which (within the framework of VA) gives a condition for finding embedded solitons, i.e., isolated 
truly localized solutions existing inside the continuous spectrum of the radiation modes. The 
criterion takes a form of orthogonality of the radiation mode in the infinite tail to the soliton core. 
To test the criterion, we have applied it to both the semi-Lagrangian truncated version of the SHG 
model and to the same model in its full form. In the former case, the criterion (combined with VA 
for the soliton proper) yields an exact solution for the embedded soliton. In the latter case, the 
criterion selects the embedded soliton with a relative error ~ 1%. 
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1 Introduction 



One of the many pioneer contributions by A.C. Newell has been the study of solitary waves (which 
we will simply call "solitons" here, in line with currently adopted terminology). These occur not 
only in integrable models, but also in many nonintegrable nonlinear- wave systems Q, §]. It is 
toward an improved understanding of solitons in nonintegrable systems, that we present this work, 
dedicated to Prof. Newell. 

In nonintegrable models, more complex forms of solitons are found as one considers higher-order 
systems. An issue of fundamental importance is to find where, in the space of the soliton parameters, 
such solutions could exist. Where this could be is dominantly determined by the model's linear 
dispersion relation, to(k). Typically, dispersion curves have gaps or partial gaps, which are values of 
to at which k would take imaginary or complex values. For example, for the nonlinear Schrodinger 
(NLS) equation, one has to = Dk 2 , where D is a positive dispersion coefficient, hence k must be 
imaginary if to < 0. This is precisely the region where one finds solitons of the NLS equation. For 
more complex systems, with several branches of the dispersion relation, one can have regions of 
to where k can be complex as well as real. Here one typically encounters "delocalized solitons" 
(see the book || and relevant examples from nonlinear-optical models in Refs. Q]), which are 
quasi-solitary waves with nonvanishing oscillatory tails. Such objects, obviously, have an infinite 
energy, and therefore are unphysical except, possibly, in finite-size systems. However, it may occur 
that the amplitude of the tail vanishes at some special values of to. Then, one has a truly localized 
object at an isolated (discrete) value of to at which a real value of k does exist. Since this value of 
to lies in a continuous part of the spectrum, these objects are called embedded solitons (ESs). A 
number of examples of ESs in physically meaningful models are now known (see examples of such 
solitons found in the hydrodynamical models in Refs. Q, and a short review in Ref. ||), their 
stability (which turns out to be semi- stability) being quite distinct from the stability of ordinary 
solitary waves 0. 

Currently, the only method for locating ESs is to search for them numerically. One does know 
that they may be found inside a partial gap of the dispersion relation, but otherwise one has no 
analytical tool for locating them (in exceptional cases, exact ES solutions can be found by guess Q). 
The first objective of this paper is to develop an approximate analytical method for locating ESs. 
The method is based on the variational approximation (VA; see a recent review of the application 
of this technique to solitons in Ref. []|] ) , which incorporates an amplitude of the infinitesimal tail of 
the related delocalized soliton as a key variational parameter. The approach will be tested on two 
versions of an ES-generating model introduced in Ref. @, which combines the second-harmonic 
generation (SHG) through quadratic [x^] nonlinearity, and the usual cubic [x®] nonlinear terms. 
The difference between the two versions of the model is that one is a full model (it has the usual 
Lagrangian representation), in terms of the x terms, while the one model is a truncated model, 
wherein, upon assuming that the fundamental-harmonic (FH) field is much stronger than the second 
harmonic (SH), some of the SH x terms are omitted. The truncated model does not admit the 
usual Lagrangian representation. However, it was this truncated version of the x' 2 ' '■ X model 
in which the above-mentioned exact analytic ES solution was found in Ref. J?]]. On the other 
hand, in the full version of the model, ES solutions could only be found by means of numerical 
methods. We will obtain a natural variational criterion which makes it possible to distinguish ESs 
from delocalized solitons in each of these systems. Furthermore, we will demonstrate that, in the 
case of the truncated model, this method yields an exact result, and in the full model, a relative 
error is ~ 1% in the prediction of the location of ES (in comparison with numerical results). 

A distinctive feature of the truncated model is that, due to a missing (omitted) term, it does not 
have a "complete" Lagrangian representation. Namely, it can be obtained from a Lagrangian, but 
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only if one term in it is not subjected to the variation with respect to the FH field, when deriving 
the system of the FH and SH equations. Accordingly, the truncated model does not conserve any 
Hamiltonian. But it is, nevertheless, a conservative system, as it conserves the norm of the solution 
(which is usually called "energy" in nonlinear optics and is different from the Hamiltonian). 

As the truncated model does not admit the full Lagrangian representation, it is necessary to 
work out a special version of the VA for it, which is another objective of this paper. In fact, this 
can be done in a very simple way: after inserting the variational ansatz into the Lagrangian and 
performing the integration in order to obtain the corresponding effective Lagrangian (the one which 
is an explicit function of variational parameters, rather than a functional depending on the field 
variables), the above-mentioned special term in the effective Lagrangian should not be varied with 
respect to variational parameters that belong to the FH component of the ansatz. Thus, other 
models similar to the above-mentioned truncated one, even though they cannot be represented in 
the usual Lagrangian form, can still be handled by means of VA. Since one cannot freely vary all 
the fields in the Lagrangian of these systems, we call them semi-Lagrangian systems. 

The rest of the paper is organized as follows. In section 2, we introduce the full and truncated 
models and the Lagrangian representations of each. In section 3, the VA for the semi-Lagrangian 
case is developed. The general VA-based analytical condition for identifying ESs (in an approximate 
form) is obtained in section 4. Section 5 concludes the paper. 

2 The full and truncated models 

Following Refs. and flffl, we first introduce the full x '■ X model: 

iu z + {l/2)u tt + u*v + 7i {\u\ 2 + 2\v\ 2 ) u = 0, (1) 

iv z - (l/2)5v tt + qv + (l/2)u 2 + 2 72 (\v\ 2 + 2|u| 2 ) v = 0, (2) 

which is written in the usual "optical" notation, so that z and t are the propagation distance and 
reduced time, u and v are the FH and SH fields, —5 is the relative dispersion coefficient at SH, q 
is the SHG mismatch, and 71 2 are Kerr coefficients. In fact, the ratio of the self-phase-modulation 
(SPM) and cross-phase-modulation (XPM) coefficients is not necessarily 1 : 2 simultaneously in 
both equations, as it is written in Eqs. (|l]) and (|2|), but this feature of the model is not a crucially 
important one. 

The form of Eqs. (p]) and (||) implies that the group-velocity dispersion is anomalous at FH, 
while at SH it may be both normal, if 5 > 0, and anomalous, if S < (both cases are physically 
possible). As for the Kerr coefficients 71 and 72, they always have one sign. Most typically, they are 
positive (corresponding to the self-focusing nonlinearity) , but may be negative too, see a detailed 



discussion of this point in Ref. [1C]. 



In many cases, the SH field is much weaker than the FH field - for instance, if the mismatch 

1 1 2 I 1 2 

is large. Then, assuming that \v\ <C \u\ , one may neglect the XPM term in comparison with the 
SPM one in Eq. (||), and the SPM term in comparison with its XPM counterpart in Eq. (|2|), which 
leads to the truncated model, 

iu z + (l/2)u u + u*v + ji\u\ 2 u = 0, (3) 
iv z - (l/2)5v tt + qv + {l/2)u 2 + ^ 2 \u\ 2 v = 0. (4) 
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In this work, we are interested in stationary fundamental-soliton solutions, which are looked in 
the form 

u(z, t) = e lkz U(t), v(z, t) = e 2ikz V(t), (5) 

where k is the FH wavenumber, and real even functions U(t) and V(t) with a single maximum at 
t = exponentially decay at t — > oo. Note that ordinary (non-embedded) solitons may exist in the 
regions 

< k < q/2, if 5 > 0; k > max{0, q/2}, if 6 < (6) 

(which implies that q must be positive for the existence of ordinary solitons if 5 is positive). On the 
other hand, ESs may exist in a range of k which does not overlap with the continuous spectrum in 
the FH equation, ([[]) or (||), but falls into the continuous spectrum of the SH equation, @ or @). 
Thus, ES may exist in the regions [cf. Eqs. (pi)] 



k > max{0, q/2}, if 5 > 0; < k < q/2, if 6 < 0. (7) 

This implies that, in the case 5 < 0, the existence of ESs makes it necessary to have positive q. 

Substituting the expressions (D into Eqs. @ and (§), we arrive at a system of ordinary 
differential equations, 

-kU + (1/2)17" + U*V + 7i (U 2 + 2V 2 ) U = 0, (8) 

-2kV - (1/2)5V" + qv + (1/2)U 2 + 2 72 (V 2 + 2U 2 ) V = 0, (9) 

where the prime stands for <i/d£. Accordingly, the system of stationary equations corresponding to 
the truncated model is 

-kU + (1/2)U" + U*V + -fiU 3 = 0, (10) 

-2W - (1/2)5V" + qV + (1/2)U 2 + A~/ 2 U 2 V = 0. (11) 
It is obvious that Eqs. (JsT) and (P) can be derived from the Lagrangian 



L = \f +0C -kU 2 - (2k + q)V 2 -\ (Uf + S - [V'f + U 2 V + ^U* + A 12 U 2 V 2 + ^V 4 

2 J —oo L 11 1 1 



dt. 
(12) 



In the case of the truncated system of stationary equations, (|l~0|) and (|Tl|), one may still obtain 
this system from the variational principle, provided we do the following. First, the last term in 
the integrand in Eq. ( |l~2| ) should be dropped. Second, the next to the last term should only be 
subjected to varying in V, but not in the U field. 

3 The variational approximation for the semi-Lagrangian system 

In this section, we focus on the application of VA to the semi-Lagrangian (truncated) system (|l0|), 



(11). The most natural variational ansatz to search for embedded solitons proper (without the tail, 



which will be considered in the next section) in this system is based on the following expressions: 

U = A sech (V2kx) x, V = B sech 2 (V2kx) , (13) 
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where the amplitudes A and B are variational parameters, while the inverse width \/2k has been 
fixed to match to the linearized form of Eqs. (§) and (H|) or (|io| ) and ( |TT1) at i — > oo. Substituting the 
ansatz ( |l~3| ) into the Lagrangian (12), dropping the last term, and then performing the integration, 
we find the effective Lagrangian of the truncated system: 



3V2kL eS = -4kA 2 



2 1 



26 



k — q 



B> + 2A 2 B + lx A 4 + - l2 A 2 B 2 . 
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(14) 



As it was said above, in order to go from the full system to its truncated counterpart, one not 
only has to omit the last term in Eq. (|i~2|), but also must take care to not vary the term 4^2U 2 V 2 
with respect to U. In terms of the effective Lagrangian (14), this means that, when deriving the 
variational equations for A and B, one should not vary the last term of the effective Lagrangian 
with respect to A. In this case, the variation with respect to A yields a simple equation which 
allows one to eliminate B, 



B = -7i A 2 + 2k. 



(15) 



Using this result, the equation produced by varying the Lagrangian with respect to B can be cast 
into the final form of a biquadratic equation for A, 



32 



7172 A 4 



64, f f 25 

l + —k l2 + 2 ll f 2 M - — 



A 2 + 4k 



2 1 



26 



0. (16) 



Thus, depending on values of the parameters, Eq. (16) may give up to two different physical 
solutions. The simplest nontrivial case (which still allows an ES to exist) is with 71 = 0, but 
72 7^ 0. Then, Eq. (|l6|) yields a single solution, 



A 



4k [2 (5 - 25) k - 5g] 
5 + 64^72 



(17) 



which is physical if it gives A 2 > 0; the amplitude B can then be obtained from Eq. (|i~5|). 

It should be stressed that, as it was demonstrated in Ref. ||], the ansatz ( |i~3| ) yields exact 
soliton solutions to Eqs. ([H]) and (11) at some uniquely selected value of k [an expression for it is 
given below in Eq. (|25|)]. Depending on the value of the mismatch parameter q, this exact solution 
may be either ES or an ordinary soliton. Comparing Eqs. ( |T5| ) and (|l^) with that solution, one 
can easily verify that expressions (|l5| ) and ( |i~6| ) are precisely parts of the exact solution. 

Although the exact solution for the ES in the truncated model is available at the single value of 
k at which it may exist, in the region @, where a continuous family of ordinary solitons is expected 
to exist, no general exact solution is known. So, to illustrate the accuracy and reliability of the 



modified VA for producing approximate solutions to the truncated equations ( |10|) and (|Tl|), in Fig. 
1 we display a typical example of a numerically found ordinary soliton in the region < k < q/2 
(in this example, 6=1), together with the analytical approximation generated by Eqs. (|l3|), (|l5|), 
and (|l6|). As one looks at the numerical solution in Fig. 1, one observes that the SH component 
goes slightly negative along the shoulder of the soliton, and then appears to oscillate as it decays. 
Note that in the region where this behavior of the SH component is observed, U is not small, and, 
according to Eq. (Q), the solution for V in this region should oscillate indeed. Thus, ordinary 
solitons in the truncated system may have fine features which are not found in the ESs. 



4 An analytical criterion to identify embedded solitons 
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4.1 General analysis 

In the case when the wavenumber k falls into the regions ([?]), ESs can exist in both the full and 
truncated models, i.e., ([[]), @ and @, (f|) 0. However, except for using the exact ES solutions 
in the truncated model to guess where ESs might exist in its full counterpart, the only known 
method for locating ESs in the full model was to search for them by means of direct numerical 
computations. 

Thus, there is a need for an analytical approach to the quest for ESs. Such an approach can 
be based on VA, if one assumes that, at values of k close to /ceSj a t which an ES exists, there also 
exists a family of delocalized solitons, with small-amplitude oscillating tails in the SH component 
that vanish when k become equal to A;es- As it immediately follows from Eqs. (^) or (11), the free 
oscillating tail of a delocalized soliton, which has an infinitesimal amplitude 6 and arbitrary phase 
shift ip, is given by the expression 

Hail = bcos (yfflS) (2k - q)t + ^) , (18) 



that neglects the nonlinear terms. Now, one may add this tail to the ansatz fll3|), to have a more 
general tailed ansatz, 

V(t) = V sol (t) + b cos (V (2k - qjt + ^ , (19) 

where V so \(t) corresponds to the ansatz for the core of the delocalized soliton. 

As b is an extra variational parameter, one should add the variational equation 

dL cS /db = (20) 

to the set of equations obtained by varying the effective Lagrangian with respect to all the other 
free parameters (irrespective of the fact if the system is complete Lagrangian or semi-Lagrangian) . 
As we are interested in the location of an ES which, by itself, has 6 = (no tail), one should set 
6 = after completing the differentiation in Eq. (|20|). This means that, prior to varying in 6, one 
should only keep terms in L c g which are linear in 6, hence Eq. ( f20|) takes the general form, 



+oo 



8L_ 
~5V 



\U=U sol (t),V=V Bol (t) 



cos V (V s ) (2k-q)t + ip )dt = 0, (21) 



with 5/5V standing for the variational derivative of the underlying Lagrangian [the one given by 
Eq. ©]. 

Let us make a couple of observations here. First, the above applies equally well to the full 
system as well as to the truncated system, since only variations in the SH field are involved, while 
the differences between the two types of systems involve solely the variation in U. Furthermore, 
since the solitons sought for are even, the expression 5L/SV with V substituted by V so \(t) is also 
even. Hence Eq. ( |2~T| ) amounts to an orthogonality condition between the infinitesimal tail and the 
soliton, 



5L\ 

J \u=u sol (t),v=v sol (t) 



cos 



(V(2/<5) (2k-q)t) dt = 0. (22) 



This consideration also shows that the phase parameter tp in the expression (|l^) is not important 
in the limit of 6 = 0. 

We notice that the above derivation circumvents the formal problem of the divergence of the 
integral expression (|D for the La grangian, when there is a tail which does not vanish as |t| — ► oo. 
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The divergence did not appear since the tail was designed as a solution to the linearized version 
of Eqs. (|) or One can readily verify that the divergence can also be eliminated in the fully 

nonlinear case, with respect to the tail (provided one suitably adjusts the frequency of the tail 
according to the nonlinearity) . 

Thus, Eq. (^) is a general criterion that can be used to locate ES solutions within the framework 
of VA. In Eq. (^), we recognize that the variational derivative (6L/SV) \u=U so i,V=V st> i(t) i s J us t the 
left-hand side (l.h.s.) of the stationary equation for SH, with U and V taken as per the chosen 
ansatz for the core of the delocalized soliton. Another point is that upon applying integration by 
parts to the second-derivative term in Eq. (11), one sees that all contributions from the linear 
terms cancel in Eq. (E3), leaving only the nonlinear terms in Eq. ( |IT| ) to determine this condition. 



4.2 Embedded soliton in the truncated system 



To test the efficiency of the criterion (22), we first apply it to the truncated system. Then, the 
variational derivative in Eq. (|22] ) is replaced by the nonlinear part of l.h.s. of Eq. (|Tl|), 



j — i 



U 2 (t) +^ 2 U 2 (t)V{t) 



cos 



(yJ(2/5) {2k-q)t} dt = 0. 



(23) 



Using the ansatz (|13|), it is easy to explicitly perform the integration in Eq. (P3j), which finally 
yields a simple result (note that the FH amplitude drops out), 



472 B 



35 -k 



2k (1 + 25) - q 



(24) 



Now, combining the above results (16) and dig), which were obtained by means of the VA, with 



the relation (24) that locates where the ES must be (and is also based on VA), one can easily verify 
that this set of three relations is precisely tantamount to the exact analytical solution for the ES 
which was found by guess in Ref. @. In particular, an eventual expression for the wavenumber of 
ES is 



k 



ES 



~(1 + 25)- 1 



q--5{A l2 + 35 7 i H 1 



(25) 



The fact that, in the case of the semi-Lagrangian truncated system, VA reproduces the exact 
ES is remarkable, although the reason for this occurrence is not fully understood. We also note 
that this exact soliton is not always an embedded one, as the wavenumber ( p5| ) is not necessarily 
restricted to the region © in which ES may exist: the wavenumber may instead fall into the region 
(|6|), in which ordinary solitons are to be found. A condition (inequality) showing whether the exact 
soliton is embedded or ordinary has already been given in Ref. [Tfl. 



4.3 Embedded soliton in the full system 



The next step is to apply the general ES-selecting criterion (22) to the full system which is based 
on Eqs. (||) an d @ ). To this end, we assume that the soliton proper may again be approximated 
by the ansatz (|T^ Then, substituting the nonlinear part of l.h.s. of Eq. (0) into Eq. ( p2f 
arrive, instead of Eq. (|24|), at a more complicated relation, 



wc 



472 B 



1 + 



B 2 2 (85 + 1) k ES ~ q 
A 2 AQ5 ■ k ES 



35 • &es 



2k ES (1 + 25) - q 



(26) 
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We test the validity of this relation in the following way: take a particular example of the ES 
in the full system that was found in a numerical form in Ref. [Q], for which 5 = 1, q = 1, and 
7i = 72 = —0.05. Borrowing values of the amplitudes A and B directly from the numerical data, 
we find that A = 3.794 and B = 2.735. Substituting these values into Eq. (f26|) yields /ces = 0.688, 
while the numerical value found in Ref. Q was &es = 0.696. Thus, the relative error of the criterion 
(H) for this case is 1.1%. 

5 Conclusion 

In this work, we have put forward two modifications to the technique of the variational approxi- 
mation for solitons. First, it may happen that a physical model does not admit the full Lagrangian 
representation, as some terms may be missing due to various reasons. In the case of the x '■ X^ 
model considered in this work, this means that one term in the Lagrangian should not be varied 
when deriving the equation for the fundamental wave. We demonstrate that the VA can be applied 
to such semi-Lagrangian systems as efficiently as to their full Lagrangian counterparts. Second, 
we have shown that, by the addition of an infinitesimal tail, which does not vanish at infinity, to 
the usual soliton ansatz, and demanding, after performing the variation, that the amplitude of the 
tail be zero, we obtain an approximate analytical criterion for locating embedded solitons inside 
a family of delocalized ones, i.e., isolated truly localized solutions existing inside the continuous 
spectrum of radiation modes. The criterion takes the form of orthogonality between the radiation 
mode contained in the infinitesimal tail and the core of the delocalized soliton. To test the criterion, 
we have applied it to both the semi-Lagrangian truncated version of the yj® : X*' 3 ' model and to 
the same model in its full form. In the former model, the criterion, combined with the VA for the 
soliton proper, yields a result which completely coincides with the previously found exact solution 
for the embedded soliton. In the latter model, the criterion predicts the wavenumber corresponding 
to the embedded soliton with a relative error ~ 1%. 
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Figure Captions 

Fig. 1. A comparison between the numerical shape (solid curves) of an ordinary (non-embedded) 
soliton of the semi-Lagrangian (truncated) system of equations (|l0|) and (|i~l"l) , as found by the shoot- 
ing method, and the result (dashed curves) obtained from the modified variational approximation 
based on Eqs. (p^), (|l5|), and (|T6|), The parameters are q = 1, S = 1, 71 = 0, and 72 = —1/4, and 
both the i ■ 1 1 1 ; ■ 1 1 ; ■ ; ' r ' 

1.2 , 
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